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ABSTRACT 

We used a combination of Hubble Space Telescope and ground based data to 
probe the dynamical state of the low mass Galactic globular cluster NGG 6101. 
We have re-derived the structural parameters of the cluster by using star counts 
and we hnd that it is about three times more extended than thought before. By 
using three different indicators, namely the radial distribution of Blue Straggler 
Stars, that of Main Sequence binaries and the luminosity (mass) function, we 
demonstrated that NGG 6101 shows no evidence of mass segregation, even in 
the innermost regions. Indeed, both the BSS and the binary radial distributions 
fully resemble that of any other cluster population. In addition the slope of 
the luminosity (mass) functions does not change with the distance, as expected 
for non relaxed stellar systems. NGG 6101 is one of the few globulars where 
the absence of mass segregation has been observed so far. This result provides 
additional support to the use of the “dynamical clock” calibrated on the radial 
distribution of the Blue Stragglers as a powerful indicator of the cluster dynamical 
age. 

Subject headings: binaries: general - blue stragglers; globular clusters: individual 
(NGG 6101) 


^ Based on observations collected at the the Very Large Telescope of the European Southern Observatory, 
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1. INTRODUCTION 

Globular Clusters (GCs) are the most populous, old and dense stellar aggregates in the 
Galaxy. They are formed by millions of stars, whose age, distance and chemical composition 
can be determined with great accuracy. For this reason GCs play a crucial role in the 
current understanding of stellar and dynamical evolution and they represent the ideal target 
to constrain the interplay between the “environment” and the evolution of stars. 

The average age of Galactic GCs (GGCs; < t >~ 12 Gyr) is typically signihcantly larger 
than the time-scales in which the internal dynamical processes occur (Meylan & Heggie 
1997). During their evolution, GCs can survive the early expansion triggered by primordial 
gas expulsion and mass loss due to stellar evolution, then their evolution is mainly driven 
by two-body relaxation thus reaching higher central concentrations and eventually the core 
collapse, while loosing stars through the boundary set by the tidal held of their host galaxy 
(see for example Heggie & Hut 2003). 

Two-body relaxation drives the long-term dynamical evolution of GCs. Because of this 
physical process, heavier objects tend to sink toward the cluster centers (mass segregation), 
while less massive stars are forced towards more external orbits. The typical time-scale in 
which two-body relaxation occurs scales with the number of stars (Spitzer 1987) and it is 
typically of the order of 1-2 Gyr in GCs (Meylan & Heggie 1997). The internal dynamics 
of stellar aggregates affect objects of any mass and it can be efficiently probed by means of 
massive test particles, like Blue Straggler Stars (BSSs), binaries and Millisecond Pulsars (e.g. 
Ferraro et ah 2001, 2003). Among them, BSSs have been successfully used for this purpose, 
since they are numerous and relatively easy to measure. Indeed, Ferraro et al. (2012) have 
shown that the BSS radial distribution can be efficiently used to rank clusters according to 
their dynamical age. Using this approach, evidence of mass segregation has been observed in 
all GCs studied so far, with only very few exceptions: u Centauri, NGC 2419, Palomar 14, 
Terzan 8 and Arp 2, where the BSS radial distribution has been found to be indistinguishable 
from that of any other population (Ferraro et al. 2006; Dalessandro et al. 2008b; Beccari et 
al. 2011; Salinas et al. 2012). In general, for clusters that haven’t reached relaxation yet, 
the radial distribution of stars of any mass is expected to be the same (apart from possible 
primordial dissimilarities). The lack of mass segregation in NGC 2419 has been conhrmed 
by Bellazzini et al. (2012) studying the the radial variation of the luminosity function (LF) 


Cerro Paranal, Chile (under proposal 091.D-0562). Also based on observations with the NASA/ESA HST 
(Prop. 10775), obtained at the Space Telescope Science Institute, which is operated by AURA, Inc., under 
NASA contract NAS5-26555. 
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of Main Sequence (MS) stars at different radio, while it has been questioned for Palomar 14 
by Frank et ah (2014) who observed an increase of the slope of the stellar mass function for 
increasing distance from the cluster center, as expected for relaxed systems. 

As a part of a large observational campaign aimed at deriving the binary fraction in 
the external regions of GGCs, we present results about the dynamical state of the low mass 
{My = —6.94; Harris 1996- 2010 edition) GG NGG 6101. This is an old {t ~ 13 Gyr, Hotter 
et ah 2010), metal poor {[Fe/H] = —1.98; Garretta et ah 2009) Galactic halo GG, with a 
low concentration (c = 0.80; Harris 1996). A recent analysis of its RRLyrae content (Gohen 
et ah 2011) revealed that NGG 6101 is an Oosterhoff type H cluster. This is consistent 
with the metallicity of the cluster, but it is unusual for its kinematical properties. In fact 
NGG 6101 is one of the very few metal-poor GGs in the Galaxy with a retrograde motion. 
Because of its peculiar kinematical properties it has been possibly connected (Martin et ah 
2004) to the Ganis Major dwarf galaxy. The stellar content of NGG 6101 has been studied 
by Sarajedini & Da Gosta (1991) and later on by Marconi et ah (2001; hereafter MOl). In 
particular MOl performed a detailed analysis of the radial distributions of different stellar 
populations of NGG 6101. They found that Horizontal Branch (HB) stars and BSS are more 
centrally concentrated than MS-Turn Off stars. They interpreted this behavior as evidence 
of mass segregation among different populations. 

In this study we investigate the dynamical state of NGG 6101 by using three diagnostics: 
the radial distribution oi (i) BSSs, (ii) MS binaries, and (in) genuine MS stars with different 
masses. According to these three indicators, and at odds with the results of MOl, we find 
that NGG 6101 does not show evidence of mass segregation. We compare this observational 
fact with theoretical expectations and the new dynamical timescale estimates obtained from 
the newly derived structural parameters. 

The paper is structured as follows: in Section 2 we present the details of our observations 
and the adopted data reduction procedures; in Section 3 we derive the main structural 
parameters of the cluster; in Section 4 the BSS radial distribution is analyzed and compared 
to that of the reference populations; in Section 5 we study the binary content of NGG 6101 
at different distances from the center and in Section 6 we analyze the radial variations of the 
luminosity and mass function of MS stars. The main results of the paper are summarized 
and discussed in Section 7. 


^The same result has been obtained by Baumgardt et al. (2009) who performed detailed comparison 
between velocity dispersion profiles and theoretical models. 
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2. OBSERVATIONS AND DATA ANALYSIS 

The data set used in this paper consists of a combination of images obtained with both 
the Hubble Space Telescope (HST) and ground based facilities (see Figure 1). 

For the analysis of the BSS radial distribution (see Section 4), we used two publicly 
available catalogs. The hrst is the HST Advanced Camera for Survey/Wide Field Camera 
(ACS/WFC) catalog published by Sarajedini et ah (2007; see also Anderson et ah 2008) in 
the context of the “ACS Survey of Galactic Globular Clusters”. This catalog samples the 
innermost ~ 120" of the cluster. We will refer to this sample as the ACS data set. The 
(V, V — I) color magnitude diagram (CMD) obtained by using the Johnson-Cousin V^round 
and /ground magnitudes, is shown in Figure 2. The second public data set used in this work 
consists of the B, V and I catalog published and fully described by MOl and obtained with 
the 1.54 mt Danish telescope at ESO/La Silla. In the following we will refer to this sample 
as the Danish sample. The corresponding CMD is shown in Figure 3. 

The binary fraction analysis (Section 5) has been performed by using deep Vhigh and 
Jbessel images obtained with the the FOcal Reducer/low dispersion Spectrograph 2 (FORS2) 
at the Very Large Telescope (Prop ID: 091.D-0562; PI: Dalessandro). The 2k x Ak pixels 
MIT Red- optimized CCD mosaic in the standard resolution mode (~ 0.25"pixel~^) was 
adopted for these observations, thus to sample the largest possible held of view (FOV; 
~ 6.8' X 6.8'). Three pointings complementary to the ACS sample have been set up (see 
Figure 1) starting from r ~ 150" (and reaching r ~ 750") from the center of the cluster. For 
each pointing, eight images in the /ressel band with exposure time Gxp = 240 sec each and 
four in Vhigh with Gxp = 510 sec, have been obtained. A dither pattern of few arcsec has 
been adopted to allow a better reconstruction of the Point Spread Function (PSF) and avoid 
CCD blemishes and artifacts. Master bias and hat-helds have been obtained by using a large 
number of calibration frames. Scientihc images have been corrected for bias and hat-held by 
using standard procedures and tasks contained in the Image Reduction and Analysis Facility 
(IRAF)|^. The photometric analysis has been performed independently for each image and 
chip (see for example Dalessandro et ah 2014) by using DAOPHOTII (Stetson 1987). For each 
frame we selected several tens of bright, not saturated and relatively isolated stars to model 
the PSF, which turned out to be well reproduced by a Mohat function (Mohat 1969). The 
parameters {a, (3) describing the PSF model have been allowed to vary with a third order 
polynomial as a function of the instrumental coordinates (x, y) within the frame. For each 
chip, the best PSF model was then applied to all sources at 2 ct above the background by using 


^Astronomy Observatory, which is operated by the Association of Universities fro Research in Astronomy, 
Inc., under the cooperative agreement with the National Science Foundation. 
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DAOPHOTII/ALLSTAR. We then created a master list of stars composed by sources detected in 
at least four frames. In the single frames, at the corresponding positions of the stars present in 
the master-list, a £t was forced with DAOPHOTII/ALLFRAME. For each star different magnitude 
estimates in each filter were homogenized and their weighted mean and standard deviation 
were finally adopted as star magnitude and photometric error (see for example Ferraro et ah 
1991, 1992). We used the stars in common with the Stetson photometric secondary standard 
catalog (Stetson 2000) to report the instrumental magnitudes to the V and I Johnson bands. 
Instrumental coordinates (x, y) have been reported to the absolute (a, S) system by using 
the stars in common with the GSC2.3 catalog and the cross-correlation tool CataXcorr0. At 
this stage the catalogs obtained for each chip are on the same photometric and astrometric 
system. They have been combined to form a single catalog that we dehned FORS2 sample. 
Stars in common between different pointings have been used to check for the presence of 
residuals in the calibration procedure. 

Both the FORS2 and the ACS catalogs have been corrected for differential reddening 
using the approach described by Milone et ah (2012). We recall the reader to this paper 
for details on the procedure. For each star, we used the 50 closest neighbors to compute 
the corresponding average differential color excess 5[E{B — V)]. We adopted extinction 
coefficients from Cardelli et ah (1989). Over the FOV covered by the three pointings, the 
total variation of E{B — V) amounts to only 0.06 mag. The resulting (V, V — I) CMDs are 
shown in Figure 2. 

Since the FORS2 data set saturates at Id ~ 18 and does not provide a complete area 
coverage (see Figure 1), we complemented it with images obtained with the Wide Field 
Imager (WFI) mounted at the MPG/ESO 2.2 m telescope to homogeneously analyze the 
density profile of NGG 6101 (Section 3). We used six long exposure images: three in the 
-Bnew band with Gxp = 120 sec (Prop ID: 069.D-0582, PI: Ortolan!) and three in V/S^esosaz 
with Gxp = 60 sec (Prop ID: 068.D-0265, PI: Ortolan!). Both the pre-reduction and the 
photometric analysis have been performed as described above for the FORS2 data set. Also 
the instrumental magnitudes have been reported to the Johnson system by using the stars 
in common with the catalog by Stetson (2000). A color equation was adopted to calibrate 
the i^NEW band, while a zeropoint was enough for the ld/89£;50843 band. Instrumental 
coordinates have been reported to the absolute system for each of the eight WFI chips using 
the stars in common with GSG2.3, as done before. The resulting (V, B — V) GMD is shown 
in the right panel of Figure 3. 


"‘CataXcorr is a code aimed at cross-correlating catalogs and finding solutions, developed by P. Montegriffo 
at INAF - Osservatorio Astronomico di Bologna, and successfully used by our group for the past 10 years. 
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We emphasize that for all catalogs, the magnitudes are in the Johnson-Cousin photo¬ 
metric system. To avoid confusion, in the following we will refer to them as “5", “W' and 

U Tit 


3. DENSITY PROFILE AND CLUSTER PARAMETERS 

We used the ACS and WF/data sets to compute the density prohle of NGC 6101 from 
direct star counts. In particular, the sample is composed of all stars in the ACS catalog 
and those in the complementary WFI data-set. We used stars with 13.5 < V < 19.5 to 
avoid incompleteness and saturation problems. As done in other works (see for example 
Dalessandro et ah 2013a), we divided the FOV in 17 concentric annuli (five of them are in 
the ACS FOV) centered on the center of gravity {Cgrav), which we have adopted to be the 
one reported by Goldsbury at ah (2010; R.A. = 16^ : 25”^ : 48®.12, Dec. = —72° : 12' : 
07".9). Each annulus has been then split in an adequate number of sub-sectors (ranging 
from two to four) according to the local density of stars and the angular coverage in the 
WFI FOV. Number counts have been calculated in each subsector and the corresponding 
densities were obtained dividing them by the sampled area. Particular attention has been 
paid to incomplete area coverage of the WFI data set starting from r ~ 1000" and to the 
inter-chip gaps. The stellar density of each annulus was then dehned as the average of the 
subsector densities, and its standard deviation was computed from the variance among the 
subsectors. We made sure to guarantee some radial overlap between the two samples in 
order to homogenize the two portions of the density prohle. The resulting surface density 
prohle is shown in Figure 4. We estimated the contribution of the Galactic held background 
by averaging the densities of the hve outermost measures, corresponding to r > 400". We 
obtain an average background density of log(pbck) ~ —2.7 stars arcsec”^. We have verihed 
that the observed background density is fully consistent with that obtained by using the 
Besancon Galaxy model simulation (Robin et al. 2003) covering an area of 1° x 1° centered 
on the position of NGG 6101 and for stars in the same magnitude limit used to build the 
density prohle. This density has been subtracted to the observed density prohle to obtain 
a “decontaminated” density distribution (solid symbols in Figure 4). As apparent from the 
hgure, the background in the external regions of NGG 6101 does not appear to be constant 
at odds with what expected, but it mildly increases as a function of the distance. This 
behavior introduces some uncertainties on the background density estimate. However we 
have checked that variations of the background density within intervals compatible with the 
observations do not ahect signihcantly the density distribution. 

We ht the radial density prohle by using an isotropic single-mass King model (King 1966) 



following the procedure fully described in Miocchi et al. (2013). The model best-htting the 
observed density prohle has a concentration c = and core radius Vc = 61. 

which yield a tidal radius r* = 1200 "l 25 o- We also derive the projected half light radius 
(effective radius) The newly determined structural parameters are in partial 

disagreement with those found in the literature and they make NGC 6101 more concentrated 
and extended than previously thought. In particular, Tc is compatible within the errors with 
the one listed by Harris (1996, Tc = 58.5"), while the values of c (and as a consequence r^) 
and Tfi are much larger than in Harris (1996; c = 0.80 and = 63"). The same result comes 
out also from the comparison with the values obtained by McLaughlin & van der Marel 
(2005), who quote c = 0.60, Tc = 59.2" and = 63.6". We argue that the difference between 
our parameters and the literature is mainly due to the fact that both Harris (1996) and 
McLaughlin & van der Marel (2005) made use of non homogeneous surface brightness prohles 
that are limited to the innermost ~ 100" from the center not allowing for an appropriate 
background subtraction. 

We have estimated the distance modulus and reddening of NGC 6101 by comparing its 
CMD to that of M 30 (NGC 7099). M 30 can be used as a template for this analysis since it 
has a similar metallicity {[Fe/H] = —2.33 ± 0.2; Carretta et al. 2009) and its distance has 
been robustly determined because of its proximity. To guarantee homogeneity, we compared 
the ACS sample of NGC 6101 with the CMD of M 30 obtained with the same instrument 
within the same survey (Sarajedini et al. 2007). 

We reported the CMD of M 30 in the absolute plane My, (H—/)o by adopting as distance 
modulus {m — M) = 14.80 and reddening E{B — V) = 0.03 (Ferraro et al. 2009). To overlap 
the CMDs and align all the main evolutionary sequences in the absolute plane (see Figure 5), 
we needed to adopt (m — M)v = 16.20 ± 0.10 and reddening E{B — V) = 0.12 ± 0.02 for 
NGC 6101, which yield to a true (unreddened) distance modulus (m — M)o = 15.83 ± 0.12 
and thus to a distance d = 14.6 ± 0.8 Kpc. We used these values throughout the paper. 

These values are slightly larger than those obtained in the literature, but in most cases 
they are still compatible within the errors. In particular, Harris (1996) report {m — M)y = 
16.07T0.1 and E{B — V) = 0.05, Sarajedini & Da Costa (1991) found {m—M)y = 16.12T0.1 
and E{B — V) = 0.06, MOl obtained {m — M)y = 16.12±0.03 E{B — V) = 0.1, while Cohen 
et al. (2011) estimated {m — M)y = 16.00±0.03 E{B — V) = 0.1 from the average luminosity 
of RRLyrae stars and ZAHB £t. 



4. THE BSS RADIAL DISTRIBUTION 


BSSs in GCs are commonly defined as those stars located along an extrapolation of 
the MS, in a region brighter and bluer than the Turn-Off (TO) point in the optical CMD 
(Sandage 1962). Their location suggests that they are more massive than the current cluster 
population (see for example Ferraro et ah 2006; Lanzoni et ah 2007a). Indeed observational 
evidence (Shara et al. 1997; Gilliland et ah 1998; Fiorentino et ah 2014) showed that BSSs 
have masses of M ~ 1.2 — 2 .OM 0 compared to a mass of ~ O.SMq for stars at the MS-TO 
of an old stellar system. 

Because of their mass, BSSs are heavily affected by dynamical friction and thus they 
can be used as natural test particles to probe the internal dynamics of stellar aggregates. 
In particular, their radial distribution has been found to probe the efficiency of dynamical 
friction (Mapelli et al. 2006; Alessandrini et al. 2014; Miocchi et al. 2015). Ferraro et 
al. (2012) suggested that the BSS radial distribution can be used to dehne the so called 
dynamical clock which is a powerful indicator of the cluster dynamical age. In this picture 
GGs with a flat BSS radial distribution are dynamically young stellar systems and they 
are defined as Family I GGs, clusters with a bimodal radial distribution have intermediate 
dynamical ages {Family IF), while those with a centrally peaked and monotonically decreasing 
BSS distribution are dynamically old {Family III). 

The BSS population of NGG 6101 has been studied by Sarajedini & Da Gosta (1991) 
and MOl. In particular MOl studied also the BSS radial distribution by using a combination 
of HST Wide Field Planetary Gamera 2 (dehned “hst” sample) and ground based data 
( “ground” sample) corresponding to the Danish sample dehned in Section 2. By using 
cumulative radial distributions and Kolmogorov-Smirnov tests, they concluded that BSSs 
have a large probability (84%) to have been extracted from the same parent population as 
HB, while HB and Tip RGB stars being more centrally concentrated than MS-TO stars. 
They interpreted these results as evidence of mass segregation. 


4.1. Population selections 

We studied the BSS population of NGG 6101 by using the ACS and Danish data sets. 
In detail, the dual sample consists of stars in the ACS catalog and those present in the 
Danish data set and complementary to the RGFFOV (see Figure 1). In this way the central 
region of NGG 6101 is homogeneously covered up to a distance r ~ 250" (corresponding 
to about 4 X Tc ) from Cgrav While the WFI data set covers a larger FOV and ensures a 
complete sampling of the cluster up to r*, we preferred this combination of catalogs for the 
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BSS analysis because of the better photometric quality of the data (see Figures 2 and 3). 

As done in previous works (see for example Ferraro et ah 1997, 2004; Lanzoni et ah 
2007a,b; Dalessandro et ah 2008a and reference therein), in order to study the BSS radial 
distribution, we need a homogeneous selection of BSSs and at least one reference stellar 
population. For the case of NGC 6101 we selected HB, RGB and MS-TO stars as reference. 

BSSs have been selected following the dehnition of stars bluer and brighter than TO. 
The selection box shown in Figure 6 follows the well defined BSS sequence brighter than 
V = 20. An additional constraint on the color ((R — I) < 0.5) has been adopted to avoid 
contamination and blends from MS and SGB stars. In this way we selected 52 BSSs, 33 in the 
ACS FOV and 19 in the complementary Danish sample. Because of the more conservative 
selection criteria adopted in this work, the BSS sample is smaller than the one obtained by 
MOl, who identihed 73 BSS adopting different magnitude and color selections. As for the 
reference populations, we selected RGB stars along the RGB and having 13.5 < V < 18.7. 
With these limits we selected a total of 510 RGB stars, 328 in the ACS sample and 182 in 
the Danish one. We selected HB stars following their well dehned sequence for (V — I) < 0.5 
to minimize the impact of held contamination. In fact, as apparent from Figure 6, the 
Galactic held describes a vertical sequence in the color range 0.7 < {V — I) < 0.9. The 
HB selection box is shown in Figure 6 and encloses 137 stars. By using the list of known 
RRlyrae stars (Filler 1981; Gohen et al. 2011; Fitzgerald et ah 2012), we identihed 11 stars 
in common, seven in the ACS and four in the Danish FOVs. The RRLyrae falling outside 
the HB selection were added to the HB sample. In this way we count 102 and 42 HB stars 
in the ACS and Danish samples respectively. In order to provide a direct comparison with 
MOl, we used also MS-TO stars as additional reference population. They have been selected 
approximately within the limits 0.65 < {V — I) < 0.80 and 20.2 < V < 20.9. It is worth 
noting that in this magnitude range, MS stars have a completeness of 100% both in the ACS 
(see Section 5.1) and in the complementary data-set (MOl). We count 3125 MS-TO stars in 
total, 2375 in the ACS and 750 in the complementary Danish FOVs. 

We obtained an estimate of the Galaxy held contamination by using the Besancon 
Galaxy model simulation (see Section 3). We counted the number of stars lying in the 
adopted (R, V — I) selection boxes. We hnd 125 and 455 Galactic held MS-TO and RGB 
stars respectively, yielding densities ~ 1 x 10“"^ stars arcsec“^ and ~ 3.5 x 10“^ 
stars arcsec”^. We don’t hnd any held star falling in the boxes adopted to select HB and 
BSS stars. 
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4.2. Radial distributions 

We first compared the BSS cumulative radial distribution to those of the reference 
populations. As done in Dalessandro et ah (2013b), we took into account the effect of 
Galactic held contamination by statistically decontaminating the RGB and MS populations. 
We divided the FOV into four concentric annuli centered on Cgrav and for each of them we 
randomly subtracted a number of stars corresponding to the average held densities quoted 
above. The decontaminated cumulative radial distributions are shown in Figure 7. As 
apparent, the four populations show very similar behaviors. Indeed when a KS test is 
applied, we obtain probabilities Rbss/rgb ~ 18%, Pbss/hb ~ 20% and Pbss/ms ~ 65% that 
they are extracted from a diherent parent population. Since they are well below the canonical 
probability limit of 95%, we can conclude that in all cases the diherences in the observed 
radial distributions are not signihcant. Also, we obtain that the radial distributions of the 
parent populations (HB, RGB and MS) are consistent with all being extracted from the 
same parent distribution. These results are in good agreement with what found by MOl 
as far as the analysis is limited to HB, RGB and BSS stars, while a diherent behavior is 
observed for MS stars. In fact, MOl found the BSS population to be signihcantly more 
centrally segregated than MS stars. In order to compare more accurately our results with 
MOl, we used also their “hst” sample and selected stellar populations using the same boxes 
dehned above. Also in this case, we do not hnd any signihcant diherence among the radial 
distribution of RGB, HB and BSS and that of MS stars. 

We also analyzed the radial distribution of the specihc frequencies Abss/Argb, Nbss/Nub 
and A^bss/^msj where A^pop is the number of stars in the corresponding population. We di¬ 
vided the FOV into four concentric annuli centered on Cgrav and in each of them we counted 
the number of BSS and that of the reference population stars. Galaxy held contamination 
has been accounted as described above. The radial distributions of the specihc frequencies 
are shown in Figure 8. In all cases, and in agreement with what observed from the cumulative 
radial distributions, we hnd that these ratios are almost constant over the entire extension 
of the ACS + Danish FOV (r < 250"). 

We computed also the double normalized ratio (RporS; Ferraro et ah 1993) for BSS 
and RGB stars. In particular, for the same annuli used before, we estimated the sampled 
luminosity using the best-ht King model discussed in Section 3, a distance modulus (m — 
M)v = 16.20 ± 0.1 and reddening E{B — V) = 0.12. We hnd that the Rpop is constant for 


^Rpop = (7Vpop/A^POp)/(l-"™P/-^tor''). where POP = RGB and BSS. is the luminosity sampled 

in each annulus and estimated by using the best-fit King model, distance modulus and reddening derived in 
Section 3, is instead the luminosity sampled in the entire field of view. 
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both RGB stars and BSSs as expected for any post-MS population (Renzini & Fusi Pecci 
1988). The same result is obtained when HB and MS stars are considered. 

Although this analysis does not cover the entire extension of NGC 6101 (it is limited to 
~ 2rfi), we note that the considered data-sets sample about the 70% of the total luminosity 
of the cluster. On the basis of what observed in tens of GGs (Ferraro et ah 2012), we do not 
expect any signihcant deviation from the general behavior at larger distances. Therefore we 
can safely conclude that the BSS radial distribution of NGG 6101 is flat (i.e., indistinguish¬ 
able from that of the reference populations). Such a distribution is not common for GGs 
and it has been observed only in few other cases so far [u Gentauri - Ferraro et ah 2006; 
NGG 2419 - Dalessandro et ah 2008b; Palomar 14 - Beccari et al. 2011; Ter 8 and Arp 2 - 
Salinas et al. 2012). As for those GGs, a flat BSS radial distribution is a strong indication 
of the absence of mass segregation in NGG 6101. 


5. THE BINARY FRACTION RADIAL DISTRIBUTION 

The fraction of binaries is an essential component of the formation and evolutionary 
processes of stellar systems. In dynamically active aggregates as GGs, binaries are thought 
to promote the formation of exotic objects like BSSs, X-ray sources and Millisecond Pulsars 
(see e.g. Me Grea et al. 1964, Paresce et al. 1992; Ferraro et al. 2006b, 2009; Heinke et al. 
2003; Pooley & Hut 2006; Xin et al. 2015). Binaries are typically more massive than the 
average single star mass in GGs (< m >~ O.SMq), therefore they tend to sink towards the 
center of GGs because of dynamical friction. As a consequence, also their radial distribution 
can be an useful tool to constrain the dynamical state of GGs. 

Up to now, three main techniques have been used to measure the binary fraction: (i) 
radial velocity variability (e.g. Mathieu & Geller 2009), (ii) search for eclipsing binaries 
(e.g. Mateo 1996; Gote et al. 1996), and (Hi) studies of the distribution of stars along the 
cluster MS in GMDs (Romani & Weinberg 1991; Bellazzini et al. 2002). In this paper we 
used the latter approach by following the method described by Bellazzini et al. (2002; see 
also Sollima et al. 2007; Dalessandro et al. 2011; Beccari et al. 2013). The basic idea 
is that the magnitude of the binary system corresponds to the luminosity of the primary 
star (more massive) increased by that of the companion. Stars on the MS obey a mass- 
luminosity relation, hence the luminosity of the binary system is a function of the mass ratio 
q = m 2 /mi of the two components (where mi and m 2 are the masses of the primary and 
secondary, respectively). Since q can assume any value between 0 and 1, in GMDs binaries 
describe a “secondary MS” in the GMD, i.e. they generate a broadening of the single stars 
MS, toward higher luminosities. 
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The binary fraction (,^) of the central regions of NGC 6101 has been estimated by 
Sollima et ah (2007) and Milone et al. (2012). In this work, we recomputed the central 
binary fraction by using the ACS sample for homogeneity, and we extended the analysis to 
the cluster peripheries by using the FORS2 data set. 


5.1. Artificial stars experiment and Galcixy contamination. 

In order to estimate the binary fraction by analyzing the “secondary MS”, a number 
of spurious effects should be considered. First, the Galactic held contamination should be 
properly taken into account. Then, a robust measures of source of broadening, like blending 
and photometric errors, should be considered. These factors are related to the quality of the 
data and can be properly studied through artihcial stars experiments. 

For the ACS sample, we used the artihcial star catalog provided with the real stars 
catalog and available at the “AGS Survey of Galactic Globular Glusters” web pag^. For 
the FORS2 sample we performed a large number of artihcial stars experiments. We followed 
the method described by Bellazzini et al. (2002; see also Dalessandro et al. 2011). We 
generated a catalog of simulated stars with an I-band input magnitude (Jin) magnitude 
extracted from a LF modeled to reproduce the observed luminosity function (LF) in that 
band and extrapolated beyond the limiting magnitude. Then to each star extracted from 
the LF we assigned a Idn magnitude by means of an interpolation along the mean ridge line 
of the cluster. 

Artihcial stars were added to real images by using the DAOPHOTII/ADDSTAR software. In 
order to avoid “artihcial crowding”, stars were placed into the images following a regular 
grid composed by 35 x 35 pixels cells (corresponding approximately to 8-9 times the usual 
FWHM of stars for these images) in which only one artihcial star for each run was allowed 
to lie. More than 250,000 stars have been simulated for the entire FOV covered by the 
FORS2 sample. The photometric reduction process used for the artihcial star experiments 
is exactly the same as described in Section 2. Those stars recovered after the photometric 
analysis have also values of Wut and lout- As already noticed in Dalessandro et al. (2011; 
see also Milone et al. 2012), the MSs of the artihcial stars GMDs are narrower than the 
observed ones. One possible reason for this ehect, is that the formal photometric errors of 
the artihcial-star catalogs systematically underestimate the true observational uncertanties. 
Another possibility is that the broadening of the observed data is due to some physical 
ehects, as the presence of multiple populations along the MS. Irrespective of the origin and 


®http://www.astro.ufl.edu/~ata/public _hstgc/ 
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as done in Dalessandro et al. (2011) we increased the formal artificial-star catalog errors 
thus to reproduce the observed error distribution as a function of magnitude. The curves of 
photometric completeness {C), dehned as the ratio between the number of stars recovered 
at the end of the procedure and the total number of stars actually simulated, are shown in 
Figure 9 for different radial bins. 

For a proper measurement of the binary fraction we performed a detailed study of the 
field contamination, using the same Besancon Galaxy model catalog described in Section 4. 


5.2. The binary fractions 

The high photometric quality and the spatial coverage provided by the combination of 
the ACS and FORS2 data sets, allow us to stu^ the binary fraction radial distribution 
from the core to up 600" from the cluster center^ We divided the FOV in four concentric 
annuli roughly corresponding to 0.0 — 1.0 r/j, 1.0 — 2.0 r/i, 2.0 — 3.5 and 3.5 — 5r/i, and 
we estimated both the minimum (.^min) and the global binary fractions (.^tot) following the 
approach described by Bellazzini et al. (2002). The analysis was performed for stars with 
19 < / < 22, where the completeness parameter iC) is safely larger than 50%. This interval 
corresponds to a mass range 0.55 < tu/Mq < 0.80 for a single star along the MS, according 
to an isochrone of 12 Gyr and Z = 0.003 (Bressan et al. 2012). 

First, we measured ^min, which is the fraction of binaries with a mass ratio gmin large 
enough to make them clearly distinguishable from single MS stars. The value of gmin depends 
on the photometric errors: we considered a color range equal to three times the photometric 
error from the MS ridge line (see Figure 3 in Dalessandro et al. 2011); in the considered 
magnitude interval and for the entire sample such a color difference corresponds to ^mm ~ 0.4. 
The contamination from blended sources and non member stars has been accounted by means 
of the artihcial star catalogs and the Besancon model simulation (Section 5.1). The radial 
distribution of is shown in Figure 10. As already observed for BSSs, binaries show 
no signihcant evidence of mass segregation. In fact, their radial distribution is almost flat 
around a mean value ^min = (6.2 ± 0.4)%. Such a value is slightly larger than what obtained 
by Milone et al. (2012) only in the central regions (gmin = 4.8 ± 0.3%). It is important to 
note, however, that Milone et al. (2012) considered only binaries with a mass ratio larger 
than ours (g > 0.5). 


^ While the FORS2 data set extends out to r ~ 800", we preferred to limit the analysis at r ~ 600" where 
the fraction of covered area coverage is more complete and the number of detected stars is appropriate for 
an accurate analysis of the binary population. 
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As a second step of the analysis, we estimated also .^tot- To this aim we created 
hundreds of synthetic CMDs including single and binary stars for different input values 
of the global binary fraction (.^in), for each run. To simulate the binaries we randomly 
extracted iVbin values of the mass of the primary from a Kroupa (2002) initial mass function 
and iVbin values of the mass of the primary component from the Fisher et ah (2005) mass 
ratio distribution, f{q), which is one of the unknowns in this kind of analysis. The impact 
of adopting different mass ratio distributions has been discussed in Sollima et ah (2007) and 
Dalessandro et ah (2011). To generate the sample of single MS stars we adopted the best £t 
to the observed mass function (MF), which is discussed in detail in Section 6. In Dalessandro 
et al. (2011) we have verihed that the adoption of different present day mass functions to 
simulate MS stars has a negligible impact on the hnal result, .^tot was then determined from 
the comparison between the artihcial and the observed CMDs: the value of providing the 
best match between the two CMDs is adopted as global binary fraction In particular, 

from the simulated catalog we computed the ratio rsim = -^biTZ-^M™ between the number 
of binary stars dehned as before and that of the synthetic MS population. The same was 
done for the observed CMD, thus obtaining rsim = -^bZnZ-^MS- 1^°^ different values of we 
computed the penalty function x^{^in) dehned as the summation of (rsim — ''"obs)^ and the 
relative probability T’(x^) was derived (Figure 11). The mean of the best-htting Gaussian 
is adopted ad best £t value of the global binary fraction .^tot- The radial distribution of 
.^tot is shown in Figure 10. Similarly to the minimum binary fraction, also .^tot shows no 
variation as a function of the distance from the center in the FOV covered in this analysis. 
This results conhrms that NGC 6101 has not experienced complete relaxation yet and does 
not show any evidence of mass segregation among its populations. 

The mean value of the global binary fraction is .^tot = (14.4 ± 0.9)%. This value 
is in good agreement with the one obtained by Sollima et al. (2007) in the ACS FOV 
(■Ctot = 15.6 ± 1.3%), while it is signihcantly larger than the one estimated by Milone et al. 
(2012; .^tot = 9-6 ± 0.6%). Note however that Milone et al. (2012) determined the value of 
.^tot assuming that it is twice the value of .^min- 


6. LUMINOSITY AND MASS FUNCTION RADIAL VARIATIONS 

To further investigate the mass segregation phenomenon in NGC 6101, we analyzed the 
radial variation of the luminosity function (LF) of MS stars. LFs give information about the 
effect of cluster internal dynamics on stars in a wide range of masses, including the faint-end 
of the MS where most of the cluster mass lies. In relaxed stellar systems, the slope of the 
LFs is expected to vary as a function of the distance from the cluster center, with indexes 
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decreasing as the distance increases, because of the differential effect of mass segregation. 

To study the LF of NGC 6101 we used the ACS and FORS2 samples. We selected 
a sample of bona-hde stars along the MS, dehned as those stars located within 3a from 
the mean ridge line, where a is the combined photometric uncertainty in the V and I 
bands. We used the same radial bins defined in Section 5 for the binary fraction analysis. 
For each radial interval we obtained the completeness-corrected LF for stars fainter than 
G = 18 and reaching the magnitude limit where the completeness parameter C is 50%. 
Completeness has been derived by means of the artificial star catalog described in Section 5.1 
(see Figure 9). The LFs obtained in the four radial bins are shown in Figure 12. We took 
the LF corresponding to the innermost radial bin as reference for comparison with the LFs 
in the other radial ranges. The line connecting the points of the reference LF is plotted as 
a solid line in Figure 12. As apparent, it well reproduces the behavior observed in the more 
external radial bins after a proper normalization. The key result of Figure 12 is that, in the 
range of magnitude and radial extension (0 < r < 5r/i) covered by our analysis, the shape of 
the different LF is clearly the same and any subtle difference lies well within the combined 
uncertainty of the LFs and of the adopted normalization. 

We derived the mass function (MF) of the same bona-hde stars used for the LF and 
within the same magnitude limits. Masses have been obtained by using the Baraffe et ah 
(1997) mass-luminosity relation. The MFs obtained in the same radial bins as the LFs are 
shown in Figure 13. The MF obtained in the ACS sample goes down to m ~ 0.15Mq, while 
it reaches at most m ~ O.35M0 in the FORS2 data set. As done before, we used the MF of 
the innermost radial bin as reference, and we normalized the more external ones to this by 
using number counts in the mass range (0.5 < tti/Mq < 0.7). As apparent from Figure 13, 
and in agreement with what obtained from the analysis of the LFs, all MFs are virtually 
indistinguishable. 


7. DISCUSSION 

Three distinct indicators, namely the radial distribution of BSSs, binaries stars and the 
LF of MS stars have been investigated to study the mass segregation in NGC 6101. We 
have traced this effect by using stars in a wide range of masses, going from ~ O.35M0, 
corresponding to the lowest mass MS stars, up to ~ 1.4M0 reached by the highest mass 
BSSs and binaries. The results clearly indicate that the cluster do not show evidence of 
mass segregation up to r ~ 600" (~ 5r/i) from the cluster center. 

Only few other GGCs are known to be in such a dynamical state: massive systems 
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(Mv < —10), namely u Centauri (Ferraro et al. 2006a) and NGC 2419 (Dalessandro et 
al. 2008b), and three GCs that are at the low mass end distribution (My ~ —5), namely 
Palomar 14 (Beccari et al. 2011), Ter 8 and Arp 2 (Salinas et ah 2012). Similarly to these 
clusters, NGG 6101 can be considered dynamically young. 

Starting from the new structural parameters obtained in Section 3, we derived the central 
(trc) and half mass radius (trh) relaxation times following the formulae by Spitzer (1987). 
We adopted the same values of {m — M)y and E{B — V) used in Section 4.2, which yield 
a distance from the Sun d = 14.6 ± 0.8 Kpc (see Section 3). We estimated the total cluster 
mass Mt by adopting an integrated V magnitude Vt = 9.2 (Harris 1996) and a mass-to-light 
ratios M/Ly = 3.17 and 2.07 (Maraston et al. 1998), appropriate for a Salpeter and Kroupa 
IMF respectively, and for a population of metallicity [Fe/H] = —2.25 and age t = 13 Gyr. 
We obtained Mt = 1.2 — 1.7 x IO^Mq which gives a central mass density pm,o = G8.5Mq 
pc“^. We thus obtained trc ~ 1-3 Gyr and Uh ~ 5.4 — 6.3 Gyr. The value of Ac is consistent 
with what reported by Harris (1996), while trh is about three times larger than the Harris 
value {trh ~ 1.7 Gyr) and more than hve times larger than the one derived by McLaughlin & 
van der Marel 2005 {trh ~ 0.9 Gyr). This is qualitatively in agreement with being larger 
than previous estimates reported in the literature. 

We note however that both trc and trh are signihcantly smaller than the age of NGG 6101 
(Gge = 13 ± 1 Gyr; Dotter et al. 2010). Therefore some evidence of mass segregation should 
be visible, at odds with the observed flat distribution of BSS and binaries, and the constant 
slope of the MF LFs. In the cases of NGG 2419 and Palomar 14, has been estimated to 
be ~ 20 Gyr (Dalessandro et al. 2008b; Sollima et al. 2011), in good agreement with the 
observational evidence of a lack of mass segregation, at least in the external regions. However 
also for NGG 2419, the central relaxation time is signihcantly smaller {trc ~ 6 Gyr) than the 
age of the cluster {t = 12 Gyr) and Ferraro et al. (2006) found for u Gentauri a relatively 
short relaxation time {trt ~ 5 Gyr), as for NGG 6101. In the case of u Gentauri a number of 
possibilities to understand this apparent disagreement between observations and theoretical 
results were examined. Gluster rotation can play a role. In fact angular momentum tends to 
keep stars out of the core, balancing the effect of mass segregation (Spurzem 2001). Another 
possibility is that u Gentauri could have been hundreds of time more massive in the past 
than observed today. Indeed it has been suggested to be the relic of a partially disrupted 
galaxy (Bekki & Freeman 2003; Tsuchiya et al. 2004). The same argument can be applied 
to NGG 2419 (van den Bergh & Mackey 2004; Mackey & van den Bergh 2005) and the 
extended tidal tails observed around Palomar 14 (Sollima et al. 2011) would suggest that it 
experienced quite recent and efficient tidal stripping events. 

For the case of NGG 6101 no evidence of internal rotation are known. However we 
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may speculate that NGC 6101 could have experienced a quite complex dynamical history. 
By means of dynamical simulations, Martin et al. (2004) suggested that NGC 6101 can be 
associated to the Canis Major dwarf galaxy accreted to the Galactic halo after an encounter 
between the dwarf and the dark matter halo of the Milky Way. Indeed the kinematical 
properties of NGC 6101 are atypical for systems with its metallicity and age. In fact it is 
one of the very few metal-poor and old GGCs with a retrograde motion (Geisler et al. 1995; 
Rutledge et al. 1997). Mackey & Gilmore (2004) suggested that the existence of old GCs 
with large cores in a relative proximity to the Galactic center (as for the case of NGC 6101) 
can be explained by the fact that they follow wide orbits in which they spend little time 
close to the Galactic center. Alternatively, they could have been accreted only recently by 
the Galactic halo. 

In this respect it is interesting to note that a linear £t to the MFs in Figure 12 limited to 
the stars with m < O.75M0 gives a power-law index a ~ —0.9. With such a value and with 
the new estimate of the central concentration obtained in Section 3 (c = 1.3), NGC 6101 
would nicely £t the correlation found by De Marchi, Paresce & Pulone (2007) according to 
which less concentrated GCs tend to have flatter MFs. This behavior was interpreted by the 
authors as an indication of a more efficient loss of low mass stars via evaporation and tidal 
stripping in less concentrated clusters. To support this possibility, evidence of tidal arms 
and distortions should be searched for this cluster. 

On the other hand, the observational facts collected in the case of NGC 2419, u Gen 
and NGC 6101 suggest that the current theoretical estimates of the central relaxation time 
are extremely rough and must be used (at least in absolute terms) with caution. On the 
other side the nice agreement among the different mass segregation indicators used in the 
case of NGC 6101 provides additional support to the use of the radial distribution of BSS 
(the so called “dynamical clock”) as a powerful indicator of the dynamical evolution of stellar 
systems (Ferraro et al. 2012). In fact, among the adopted mass-segregation indicators, BSSs 
are significantly brighter than the others (MS stars and binaries), their analysis is simpler 
and less prone to observational bias (as completeness) and assumptions. Hence the BSS 
radial distribution represents the clearest indicator of mass segregation in stellar systems. 


This research is part of the project COSMIC-LAB (http://www.cosmic-lab.eu) funded 
by the European Research Council (under contract ERC-2010-AdG-267675). 
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Fig. 1.— Schematic map of the entire database used in this work with respect to the position 
of the center of the cluster. The solid and dashed circles represent the clusetr core and tidal 
radii respectively. 
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Fig. 3.— Left panel: (V, V — I) CMD of the Danish data set published by Marconi et ah 
(2001). Here only the stars lying ontside the ACS FOV are shown. Right panel: (H, B — V) 
CMD obtained from the WFI sample. In order to minimise the contamination from Galactic 
held stars, we show here only stars located at r < 500" from the cluster center. 
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Fig. 4.— Observed star count density profile as a function of radius (open squares). The 
dashed line represent the density value of the Galactic held background, obtained averaging 
the hve outermost points. The black hlled dots are densities obtained after background 
subtraction (Secion 3). The best-ht single-mass King model is also overplotted to the ob¬ 
servations (solid line). The structural parameters are labeled. The lower panel shows the 
residuals between the observations and the best-ht model. 
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Fig. 5.— Superposition of the CMD of NGC 6101 (black open circles) to that of M 30 (red 
crosses) in the absolute My, (V — I)o plane. For M 30 we adopted the values of distance and 
extinction obtained by Ferraro et ah (2009). A distance modulus (m — M)y = 16.20 ± 0.1 
and a color excess E{B — V) = 0.12 were then adopted for NGC 6101 to match the main 
evolutionary sequences of M 30. 














Fig. 8.— From top to bottom, radial distributions of the double normalized ratios of the 
BSS (i?Bss; solid circles) and RGB stars (Rrgb; grey regions), and the specihc frequencies 
Nbss/Nrgb, Nbss/Nub and Nbss/Nms specihc frequencies (solid circles). 





Fig. 9.— Photometric completeness C as a function of the I band magnitude for the ACS 
and the FORS2 data sets and for different radial bins. 
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Fig. 10.— Minimum and global binary fractions as a function of the radial distance from 
the cluster center. The dashed lines represent the mean value, while the dotted ones dehne 
the l(j level. 



















Fig. 12.— LFs in the / band obtained from the ACS and the FORS2 data sets in different 
radial bins. Measurements are shifted by an arbitrary amount to make the plot more read¬ 
able. The solid line is the LF obtained obtained in the ACS sample for r < r^. It has been 
overplotted to the outer LFs for comparison. 
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Fig. 13.— MFs derived by using the Baraffe et al. (1997) mass-luminosity relation. Radial 
bins and symbols are the same as in Figure 12. 





